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Abstract 


We study in detail the phase structure of a holographic p-wave superconductor 
model in a five dimensional Einstein-Maxwell-complex vector field theory with a 
negative cosmological constant. To construct complete phase diagrams of the model, 
we consider both the soliton and black hole backgrounds. In both two cases, there 
exist second order, first order and zeroth order phase transitions, and the so-called 
“retrograde condensation” also happens. In particular, in the soliton case with the 
mass of the vector field being beyond a certain critical value, we find a series of phase 
transitions happen such as “insulator /superconductor/insulator/superconductor” , as 
the chemical potential continuously increases. We construct complete phase diagrams 
in terms of temperature and chemical potential and find some new phase boundaries. 
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1 Introduction 


The strong/weak duality characteristic of the anti-de Sitter/conformal field theory cor- 
respondence (AdS/CFT) [1, 2, 3] enables us to study the properties of strong coupled 
systems by a weak coupled AdS gravity with one higher dimension. There are a lot of 
strong interacting phenomena in condensed matter physics, which are thought to be a 
good place where the AdS/CFT correspondence is applicable (see refs. [4, 5, 6, 7] for re- 
views). In order to describe more general phenomena in this framework, recent efforts 
are devoted to generalizing the correspondence to those systems with less symmetries (see 
refs. [8, 9, 10, 11, 12, 13, 14] for example) and describing the far-from equilibrium problems 
quantitatively (see refs. [15, 16, 17, 18, 19] for example). 

One of the most studied objects is the holographic superconductor (superfluid). In the 
simplest model [20, 21], the onset of superconductivity is characterized by the instability 
to form complex scalar hair in the bulk black hole spacetime, which corresponds to conden- 
sation of composite charged scalar operator spontaneously breaking U(1) symmetry below 


a critical temperature in the CFT side. This toy model describes a holographic s-wave su- 
perconductor/conductor phase transition, i.e., a transition from black hole with no scalar 
hair (normal phase/conductor) to the case with scalar hair at low temperatures (supercon- 
ducting phase). By using AdS soliton as the background in this Einstein-Maxwell-complex 
scalar theory, the authors of ref. [22] constructed a holographic model mimicking a s-wave 
superconductor/insulator transition at zero temperature. The normal insulating phase 
is described by a pure AdS soliton which exhibits a mass gap and is used to describe a 
confining phase in a dual theory [23]. Adding a chemical potential u to the bulk, the 
AdS soliton will become unstable to developing scalar hair for sufficiently large chemical 
potential, which describes a superconducting phase. It was shown that as one changes 
the temperature T and the chemical potential u there are as many as four phases in this 
model, including the AdS soliton, the AdS black hole and their superconducting phases. 
The complete phase diagram in terms of T and u was constructed in ref. [24]. 

Recently, the authors of ref. [25] proposed a holographic p-wave superconductor model 
by introducing a complex vector field p, charged under a Maxwell field A, in the bulk, 
which is dual to a strong coupled system involving a charged vector operator with a global 
U(1) symmetry. This setup meets the minimal requirement to construct a holographic 
p-wave superconductor model, since the condensate of the dual vector operator may spon- 
taneously break the U(1) symmetry as well as spatial rotation symmetry. Actually, worked 
in the probe limit neglecting the back reaction of the matter fields, the vector hair appears 
at low temperatures and the condensed phase exhibits an infinite DC conductivity and a 
gap in the AC conductivity. In this sense, this Einstein-Maxwell-complex vector model can 
be regarded as a holographic p-wave model. t The probe approximation is only justified 
in the limit of large q with gp, and qA, fixed. In the paper [28], the full back reaction 
was taken into account in the four dimensional black hole background. Depending on the 
charge g and mass square m? of the complex vector field, the model presents a rich phase 
structure. One can find zeroth order, first order and second order phase transitions. There 
is also a “retrograde condensation” in which the hairy black hole solution exists only for 
temperature above a critical value and is thermodynamically subdominant. In particu- 
lar, we can see in this model the conductor/superconductor/conductor phase transitions 
when the temperature is continuously lowered. In addition, the behavior of entanglement 
entropy in this model was discussed in ref. [29]. 

The Einstein-Maxwell-complex vector model can be directly generalized to the AdS 
soliton background to describe a holographic p-wave superconductor /insulator phase tran- 
sition. Actually, in the probe approximation, the vector condensate and instability induced 
by a magnetic field in a (4+1) dimensional AdS soliton background was studied in ref. [30], 
where a vortex lattice structure can form in the spatial directions perpendicular to the 
applied magnetic field.? Note that in ref. [30] the vector condensation and AdS soliton 


'The first holographic p-wave model was given in ref. [26] by introducing a SU(2) Yang-Mills field into 
the bulk, where a gauge boson generated by one SU(2) generator is dual to the vector order parameter. An 
alternative holographic realization of p-wave superconductivity comes from the condensation of a 2-form 
field in the bulk [27]. 

?In the ref.[31], a similar phenomenon that the magnetic field leads to the instability was also found in 


instability can be totally induced by an applied magnetic field because of the existence of 
the coupling between the magnetic moment of the vector field and the background mag- 
netic field in this model. Furthermore, it was found that this toy model is a generalization 
of the SU(2) p-wave model [26] in the sense that the vector field has a general mass and 
gyromagnetic ratio. 

As we have shown in the black hole case [28], it is worthwhile to consider the full 
back reaction of the matter fields on the soliton geometry and to find all possible phase 
behaviors. This is just one of the purposes of the present paper. On the other hand, 
in order to construct complete phase diagrams of the model, we will also study phase 
behaviors of the model in a (4+1) dimensional black hole background by generalizing the 
study in ref. [28]. Since in this paper we do not turn on magnetic field, the model is 
left with two parameters, i.e., the charge q and mass square m? of the complex vector 
field p,. More precisely, m? determines the scaling dimension of the dual vector operator, 
while q controls the strength of the back reaction. Depending on q and m?, the model in 
the black hole case exhibits all known phase structure reported in ref. [28]. In the soliton 
geometry, we can find second order transition, first order transition, zeroth order transition 
as well as “retrograde condensation” in the sense that the hairy soliton solution appears 
only for chemical potential below a critical value and has free energy larger than the pure 
AdS soliton. Those phase behaviors can be comparable to those in the black hole case. 
Nevertheless, there is an additional interesting behavior in the soliton case. In some region 
of model parameters, the condensate will be absent in some range of chemical potential, 
where the condensate would be expected to appear naively. This leads to a series of phase 
transitions. In a typical example, as one increases the chemical potential, we will encounter 
for the insulator /superconductor /insulator/superconductor phase transitions. 

Taking both the soliton geometry and black hole case into account, we obtain four 
distinct solutions: AdS soliton, AdS black hole and their corresponding hairy solutions. 
We construct the full phase diagram as a function of chemical potential and temperature for 
various q and m?. Due to the rich phase behaviors of the model, the T-y phase diagram is 
much more complicated and interesting than the holographic s-wave model [24] and SU(2) 
p-wave model [33, 34]. We can find some new phase boundaries that are never reported in 
the literature in the framework of AdS/CFT correspondence. 

The organization of this paper is as follows. In the next section, we introduce the 
holographic model and give no-hair soliton and black hole solutions. In section 3, we 
give our ansatz for the hairy soliton solution and specify the boundary conditions to be 
satisfied. Then the phase transition for various parameters will be discussed in detail. In 
section 4, a parallel discussion for the case of black hole geometry will be presented. The 
complete phase diagrams in terms of temperature and chemical potential are constructed in 
section 5. The conclusion and some discussions are included in section 6. Some calculation 
details are shown in appendixes. 


a SU(2) model in the AdS Schwarzschild background and hard wall cutoff model. 


2 The Holographic Model and No-hair Solutions 


Let us begin with a (4+ 1) dimensional Einstein-Maxwell-complex vector field theory with 
a negative cosmological constant [25] 


a= aa | ËV IR+ = + Lm), 
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where L is the AdS radius which will be set to be unity and K? = 87G is related to the 
gravitational constant in the bulk. g is the determinant of the bulk metric gv and p, is 
the complex vector field with mass m and charge q. We define Fiw = Vp Ar — V, A, and 
Puv = Dapo — Dopp with the covariant derivative D, = V, —iqA,. The last non-minimal 
coupling term characterizes the magnetic moment of the vector field p,, which is crucial 
in the presence of an applied magnetic field [25, 30]. However, in the present study, since 
we only consider the case without external magnetic field, this term will not play any role. 
So this model is left with two parameters m? determining the scaling dimension of the 
dual operator and q controlling the strength of the back reaction of matter fields on the 
background geometry. 
The full equations of motion deduced from the action (1) are the complex vector field 
equations 
D” py — py + iqyp” Fin = O, (2) 


Maxwell’s equations 
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as well as the gravitational field equations 
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Since p,, is charged under the U(1) gauge field, according to AdS/CFT correspondence, its 
dual operator will carry the same charge under this symmetry and a vacuum expectation 
value of this operator will then trigger the U(1) symmetry breaking spontaneously. Thus, 
the condensate of the dual vector operator will break the U(1) symmetry as well as the 
spatial rotation symmetry since the condensate will pick out one direction as special. 
Therefore, viewing this vector field as an order parameter, the holographic model can 
be used to mimic a p-wave superconductor (superfluid) phase transition. The gravity 
background without vector hair (p,, = 0)/with vector hair (p, # 0) is used to mimic the 
normal phase/superconducting phase in the dual system. 


First, let us consider solutions with p, = 0 and A; = ¢(r). One solution with planar 
symmetry is the AdS Reissner-Nérdstrom (RN-AdS) black hole, which reads [4] 
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The horizon locates at r, and pu is the chemical potential of the dual field theory. The 
temperature T of the black hole is 


2 
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In the following, we will assume the coordinate z is compactified with period I’, whose 
value will be given below, while x and y are in (—oo, 00). There is another trivial planar 
solution, the so-called AdS soliton, which is given by 


ds? = —f (r)dt? + + r?(dx? + dy? + dz’), 


(5) 
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where ro is the tip of the soliton. To avoid the potential conical singularity at r = rọ, the 
spatial direction 7 must make an identification with period 


r= (8) 


T 
>, 
Thus, the AdS soliton is cigar shaped with the asymptotical geometry Rt? x St near the 
AdS boundary r — oo. Because the spacetime exists only for r > rọ, the dual field theory 
is in a confined phase and has a mass gap, Eg ~ ro. Since the time component of metric 
is regular at the tip, the soliton can be associated with any temperature in the Euclidean 
sector. 

In the holographic setup, the RN-AdS black hole (5) corresponds to a conductor while 
the AdS soliton (7) describes an insulator [22]. In order to obtain a superconducting phase 
in both soliton and black hole backgrounds, we need to find solutions with non-trivial p, 
in the bulk. That is what we will do in the following. 


3 AdS Soliton with Vector Hair 


3.1 Equations of motion and boundary conditions 


To construct homogeneous charged solutions with vector hair in the soliton background, 
we adopt the following ansatz 


= —— + r?2(—f(r)dt? + h(r)dz? + dy? + g(r)e-X dn), F 


pdz” = pz(r)ìdz, A, dz” = ọ(r)dt, 


where g(r) vanishes at the tip r = ro of the soliton. Further, in order to obtain a smooth 
geometry at the tip ro, 7 should be made with an identification 
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This gives a dual picture of the boundary theory with a mass gap, which is reminiscent of 
an insulating phase. 

One finds that the r component of (3) implies that the phase of py must be a constant. 
Without loss of generality, we can take py to be real. Then, the independent equations of 
motion in terms of the above ansatz are deduced as follows 
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where the prime denotes the derivative with respect to r. 

The full coupled equations of motion do not admit an analytical solution with non-trivial 
Pz- Therefore, we have to solve them numerically. We will use shooting method to solve 
equations (11). In order to find the solutions for all the six functions F = { pz, Q, f, g, h, x} 
one must impose suitable boundary conditions at both conformal boundary r — oo and 
the tip r = ro. 


In order to match the asymptotical AdS boundary, we have the general falloff near the 
boundary r — oo as 
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where the dots stand for the higher order terms in the expansion of 1/r and Ay = 1 
V¥1+m?. 3 In general, in the above expansion we must impose pr_ = 0, which meets 
the requirement that the condensate appears spontaneously. According to the AdS/CFT 
dictionary, up to a normalization, the coefficients u, p, and pz, are regarded as chemical 
potential, charge density and the x component of the vacuum expectation of the vector 
operator J” in the dual field theory, respectively. 

We impose the regularity conditions at the tip r = rg, which mean that all our functions 
have finite values and admit a series expansion in terms of (r — ro) as * 


F= F(ra) + F (ra) -rot (13) 


By plugging the expansion (13) into (11), one can find that there are six independent 
parameters at the tip {ro, Pz(ro), (ro), f (ro), A(ro), X(ro)}. However, there exist four useful 
scaling symmetries in the equations of motion, which read 


X>xX+A, nen, (14) 
$> àd, t>, Fors, (15) 
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where in each case À is a real positive constant. 

By using above four scaling symmetries, we can first set {ro = 1, f(ro) = 1, h(ro) = 
1, x(ro) = 0} for performing numerics. After solving the coupled differential equations, one 
should use the first three symmetries again to satisfy the asymptotic conditions f(oo) = 1, 
h(oo) = 1 and y(oo) = 0. 5 We choose (ro) as the shooting parameter to match the 


3The m? has a lower bound as m? = —1 with A, = A_ = 1. In that case, there exists a logarithmic 
term in the asymptotical expansion of px. One has to treat such a term as the source set to be zero to 
avoid the instability induced by this term [37]. In this paper, however, we always consider the case with 
m? > -—1. 

4In some extremal cases, the functions do not admit such Taylor series expansion. In appendix B, it is 
shown that a type of solutions exists with a finite nonzero I’, but is expressed in terms of fractional order of 
(r—ro). This type of solutions corresponds to the limit case of the discontinuous points in subsection 3.3.1 

5Here we assume the difference of y(r) between at the tip and at the boundary is finite. However, it 
does not hold in the limit case of the discontinuous points discussed in subsection 3.3.1. The details will 
be discussed in appendix B 


source free condition, i.e., py_ = 0. Finally, for fixed m? and q, we have a one-parameter 
family of solutions labeled by p,(ro). After solving the set of equations, we can read off 
the condensate (Jey, chemical potential u and charge density p from the corresponding 
coefficients in (12). It should be noticed that different solutions obtained in this way will 
have different periods I for 7 direction. We should use the last scaling symmetry to set all 
of the periods Į equal in order to obtain same boundary geometry. We shall fix I to be m 
in this paper. 


3.2 Free energy 


In the soliton background, there are two kinds of solutions: one is the pure AdS soliton 
without the vector hair and the other is the hairy one. To determine which kind of solutions 
is thermodynamically favored, we must calculate the grand potential Q of the system in the 
grand canonical ensemble. In gauge/gravity duality the grand potential Q of the boundary 
thermal state is identified with temperature times the on-shell bulk action in Euclidean 
signature. One should also consider the Gibbons-Hawking boundary term for a well-defined 
Dirichlet variational principle and further a surface counterterm for removing divergence. 
Since we consider a stationary problem, the Euclidean action is related to the Minkowski 
one by a minus sign as 


— 2K? S Fuclidean = | Pavar +12 + Lm) + f drv -h(2K — 6), (18) 
Too 
where g is the determinant of the bulk metric (one should not confuse it with the function 
g(r) appearing in the metric ansatz), h is the determinant of the induced metric on the 
boundary, and K is the trace of the extrinsic curvature. © By using the equations of 
motion (11) and the asymptotical expansions of matter and metric functions near the 
boundary, the grand potential Q turns out to be 


2K NQ 
V3 


= ga + 4( fa + ha — Xa), (19) 


where V3 = f dxdydn. There is a little subtle issue dealing with soliton case, since the soli- 
ton background has no horizon and the associated Hawking temperature vanishes. Nev- 
ertheless, one can introduce an arbitrary inverse temperature 1/T as the period of the 
Euclidean time coordinate. The integration over the Euclidean time in the Euclidean ac- 
tion just gives the factor 1/T, which cancels the temperature factor in the Euclidean action 
and leads to a finite grand potential. In the case of m? > —1, it is found that f4+h4— y4 = 0 
(see appendix C). Since we have scaled I to be 7, g4 = —1 for the pure AdS soliton solution, 
namely, in the normal insulating phase. 


In principle, we should also consider the surface counterterm for the charged vector field p,, but one 
can easily see that this term makes no contribution under the source free condition, i.e., pr- = 0. 


3.3 Phase transition in AdS soliton backgrounds 


In what follows we will look for hairy soliton solutions numerically. Our strategy is to fix 
m? and for each m? we scan a wide range of q which determines the strength of the back 
reaction of matter fields on the background. The precision of numerical calculation limits 
us to investigate the entire range of q. Nevertheless, our numerical calculation reveals 
that the system exhibits qualitatively different behavior depending on concrete value of 
m*. More specially, there are two critical mass squares, m2, = 0.218 + 0.001 and m2. 
Although it is difficult to determine the exact value of me numerically, it is suggested to 
be zero by numerical analysis. Depending on q and m?, we can find second order, first 
order and zeroth order phase transitions. There also exists the “retrograde condensation” 
in which the hairy soliton solution exists only for chemical potential below a critical value 
and is thermodynamically subdominant. What’s more, a discontinuous condensed phase 
appears when m? > m2, and q is less than some critical value, which is a new kind of 
phase transition that does not exist in other models. We shall give more details in the 
following. Since we would meet with more than one critical chemical potential at which a 
phase transition happens as we increase the chemical potential, for brevity, we shall denote 
the first critical chemical potential as e1, the second critical chemical potential as 1.2, and 
so on. 


3.3.1 Phase transition for m? > m2, 


For the case m? > m3, there exist three critical charges qı, q2 and q3 with q1 < q2 < q3. 


Their values depend on m? we choose. The three critical charges divide the parameter 
space of q into four regions, shown in table 1. The phase behavior changes qualitatively 
in each region. For the case q > q3, the condensed phase will appear above ue through 
a second order transition and becomes thermodynamically preferred. As we decrease q to 
a value smaller than q3, the order of transition from the normal (insulator) phase to the 
superconducting phase becomes first order. Much more interesting thing happens when 
we continue decreasing the value of q past q2. Starting from the small u region, the system 
undergos a first order transition from the normal phase to the condensed phase at uc, and 
as u increases to Hep there is a zeroth order transition back to the normal phase, then at 
the larger chemical potential ue3, it comes back to the condensed phase again by a zeroth 
order transition. There is a region [leg < H < Heg in which hairy solutions do not exist, so 
the condensed phase should jump to the normal phase. As one can see clearly, there are 
three phase transitions as we increase the chemical potential. The later two zeroth order 
phase transitions have not been reported in any holographic model. For the case q < q, 
the system undergoes a zeroth order transition from the normal phase to the condensed 
phase at u«a. The phase transition and their orders are summarized in table 1. 

As a typical example, we consider the case m? = 5/4. We find similar results for other 
values of m? in the region m? > m?2,. For m? = 5/4, the three critical values of charge are 
q3 ~ 1.5600, go ~ 1.5345, and qı œ 1.5020, respectively. The pure AdS soliton survives for 
arbitrary value of q. Nevertheless, we can find additional solutions with non-vanishing px, 
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Charge Phase transition and its order, m* >m 
| gas | wiper, S | per, 2 | > ber, SC 


Hcl < H < Hez, SC 
Table 1: The phase transition and its order with respect to the charge and chemical 
potential in the case of m? > m2,. In the table, S=normal phase/insulator with pure AdS 
soliton, SC=condensed phase/superconductor with hairy soliton solutions. 0%", 1% and 2°4 
stand for the zeroth order, first order and second order phase transitions, respectively. 


which are thermodynamically preferred for sufficiently larger chemical potential. That is 
to say, for each value of g, there must be a phase transition from the normal phase to the 
condensed phase occurring at a certain chemical potential. From the perspective of dual 
system, it means that a charged vector operator obtains a non-zero vacuum expectation 
value (J,) # 0, which breaks the U(1) symmetry as well as the rotation symmetry in x — y 
plane spontaneously. ” 

For q > q3, we take q = 1.600 as a typical example. Figure 1 shows the grand potential 
and condensate as a function of chemical potential, from which one can see that (.J,) arises 
continuously from zero at uea. The grand potential Q in the left plot of figure 1 indicates 
that above ue the configuration with non-vanishing vector “hair” is indeed thermodynam- 
ically preferred to the normal phase. It is a second order phase transition with the critical 
behavior (J,) ~ (u/He1 — 1)'/? near the critical point. 

For the case qo < q < qs, such as q = 1.540 in figure 2, the transition from the 
normal phase to the condensed phase becomes first order. As one can see in figure 2, the 
condensate with respect to chemical potential is multi-valued and the free energy develops 
a characteristic “swallow tail”. The condensate has a jump from zero to the upper branch 
of the hairy solution at uea œ 1.8054. 

There is an additional complication when we decrease the value of q past q2. A concrete 
example for g = 1.534 is presented in figure 3. As one increases the chemical potential, 
there is a first order transition at uea ~ 1.8094. However, if one continues increasing the 
chemical potential, there will be two discontinuities at 2 ~ 1.9085 and ue3 ~ 1.9499, 
between which there does not exist any hairy solution. So there is a kind of zeroth order 
phase transition from the condensed phase to the normal phase at {4.2 and a kind of zeroth 
order phase transition from the normal phase to the condensed phase at pu-3. These kinds 
of zeroth order phase transitions do not come from thermodynamics but from the kinetics 
of the field equations (11). The details will be shown in appendixes A and B. 

Recall that the pure AdS soliton resembles an insulating phase while the hairy soli- 
ton mimics a superconducting phase, our result suggests that there is a kind of “insu- 
lator /superconductor/insulator/superconductor” phase transitions by increasing chemical 


7The breaking of the rotation symmetry in x — y plane is due to the fact that (J,) chooses x direction 
as special. However, this anisotropy does not display in the stress energy tensor of the dual field theory. 
In fact, the xx-component and yy-component of the stress energy tensor are equal. For more details one 
can see in appendix C. 
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Figure 1: The grand potential and vacuum expectation value of vector operator with 
respect to the chemical potential in the case of m? = 5, q = 1.600. In the left plot, the 
grand potentials of the pure soliton and hairy soliton are described by the green dashed 
line and blue solid line, respectively. The condensate begins at uea œ 1.7070. The right 
plot shows the condensate of the vector operator (J,) 4 0 when u > pa. 
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Figure 2: The grand potential and condensate with respect to the chemical potential in 
the case of m? = 5, and q = 1.540. In the left plot, the grand potentials of the pure 
AdS soliton and hairy soliton are denoted by the green dashed line and blue solid line, 
respectively. The condensate begins at ue ~ 1.8045. In the right figure, the blue solid line 
stands for the physical part of the condensate value for the condensed phase, which shows 


the condensate has a jump from zero to the upper branch when p > pa. 


12 


1k 
—— a — 
osk i : 
<j 927 
x 0.6 i 
0.4- i 
0.2 ; 
f3 1.85 1.9 1.95 2 
u 


Figure 3: The grand potential and condensate with respect to the chemical potential in the 
case of m? = 5 and q = 1.534. In the left plot, the grand potentials of the pure AdS soliton 
and hairy soliton are denoted by the green dashed line and blue solid line, respectively. 
The condensate begins at uce ~ 1.8094. However, there does’t exist the condensed solution 
when the chemical potential is in between uez ~ 1.9085 and e3 ~ 1.9499 (The red solid 
line region in the figure). The right plot shows the condensate has a jump from zero to the 
upper branch at u = uea. When the chemical potential is in the region with the red solid 
line, the condensate falls into zero again. 


potential. It is interesting to see whether there is any realistic material that can display 
such a phase structure. 

As we continuously decrease q, the value of ue decreases while the value of uez increases. 
There exists a critical value of charge qı ~ 1.5020 at which the first order phase transition 
point ue coincides with the first zeroth order point e2. In this case, the first part of 
condensed phase such as in figure 3 disappears (see figures 4 and 5). When the value of q 
is less than q4, the first order transition point disappears and the first discontinuous point 


Leg appears above the line 2K7Q/V3 = —1, thus is un-physical. This case for q < q is 
shown in figure 4 with q = 1.490. There is only a zeroth order transition at critical point 
Ley ~ 1.9499. 


The (u-q) phase diagram for the case with m? = 5/4 is drawn in figure 6. One should 
note that the horizontal axis and vertical axis in the figure are labeled by 1/q and qu, 
respectively. It is convenience to compare with the results in the probe limit, where we 
take the limit q —> oo keeping qp,, and qA, fixed. Therefore, in order to compare our results 
to the case in the probe analysis, we should make the scaling transformation pp — py, 
and A, — qA,. Under such transformation, the chemical potential u becomes qu and 
the temperature T becomes T/q. Therefore, they are the quantities qu and T/q that are 
comparable to those in the probe limit. One can see from figure 6 that for sufficiently large 
q, the critical value of qu approaches to a constant, which quantitatively agrees with the 
probe analysis in appendix D. Furthermore, let us mention that the zeroth order phase 
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Figure 4: The grand potentials with respect to the chemical potential in the case of m? = 
5, q = qı œ~ 1.502 (left plot) and q = 1.490 < qı (right plot). The grand potentials of the 
pure AdS soliton and hairy soliton are denoted by the green dashed line and blue solid 
line, respectively. In the case with q ~ 1.502, the first point of discontinuity of the grand 
potential coincides with the first order phase transition point. The right plot shows the 
case of q < qı. Note that only the part which has the lowest free energy is physical. 
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Figure 5: The condensate of the vector operator in the case of m? = 5, q = qı œ 1.502 
(left plot) and q = 1.490 < qı (right plot). For the condensed phase, only the blue line 
is thermodynamically favored, while the dashed orange part is un-physical. Thus, the 
condensate has a sudden jump form zero to a nonzero value. There exists a zeroth order 
phase transition from the normal phase to the condensed phase. 
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Phase diagram (m?=5/4) 
32 T T r 
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3r Condensed Phase 


Normal Phase 


2.7F Normal Phase 


Figure 6: The phase diagram for the case with m? = 5/4. The red line stands for the 
second order phase transition, the green line stands for the first order phase transition, 
and the black line stands for the zeroth order phase transition. Numerical results show 
uq œ 2.7852 when 1/q is close to zero. 


Phase transition and its order, m2, < m? < må 
p> Her, SC 


a> fa SC 
i> fia SC 


Table 2: The phase transition and its order with respect to the charge and chemical 
potential in the case of m2, < m? < m?2,. In the table, S=normal phase/insulator with pure 
AdS soliton solution, SC=condensed phase/superconductor with hairy soliton solution. 
0 1%, and 24 stand for the zeroth order, first order and second order phase transition, 


respectively. 


transition can only appear in the strong back reaction case. Finally in some range of q, 
one can first see one first order transition and then two zeroth order transitions with the 
increase of chemical potential. 


3.3.2 Phase transition for m2, < m° < må 


In the case of m? 2 , we find two critical values of q denoted as qa and qe with qa < 


qə. For small back reaction with q > q», there is a second order insulator/superconductor 
transition at some critical chemical potential j1.;. As we decrease q to the range qa < q < Q, 
the order of transition becomes first order. For larger strength of back reaction with q < qa, 
the condensed phase appears through a zeroth order phase transition. The main results 
are shown in table 2. The concrete values of qa and q depend on m?. To give more details, 
we shall consider a typical example with m? = 0 in this subsection. The two critical values 
of q for m? = 0 is qa ~ 1.172 and q œ 1.317. 

For q > q, the situation is very similar to the case with q > q3 in m? = 5/4. There 


o LM <m 
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Figure 7: The grand potentials with respect to the chemical potential in the case of m? = 0, 
q = 1.500 (left) and g = 1.200 (right). In both figures, the grand potentials of the pure 
AdS soliton and hairy soliton are described by the green dashed line and blue solid line, 
respectively. In the left plot, the condensate begins at uea ~ 1.5103, while in the right plot 
the condensate begins at ue ~ 1.8207. 


exists a critical chemical potential ue at which a second order transition happens. As we 
go on decreasing q past q, the transition from the normal phase to the condensed phase 
is first order. The behaviors of the grand potential and condensate are shown in figures 7 
and 8, respectively. 

When q is slightly below q, the curve of condensate (or grand potential) is continuous, 
while as q is a little larger than qa, the curve becomes discontinuous. In the previous 
subsection, this discontinuity leads to two zeroth order transitions for qi < q < q (see 
figure 2) and one zeroth order transition for q < q (see figure 4). In the present case with 
q => qa, although there exist two discontinuous points in the grant potential curve denoted 
as upper point with larger grand potential and lower point with smaller grand potential, 
as shown in the left plot of figure 9, both points appear in the region where the grand 
potential of the condensed solution is larger than the one in the normal phase, thus are 
thermodynamically disfavored and will not appear in the physical phase space. So we are 
left with only a first order transition. The lower point moves down as q decreases and 
finally its grand potential becomes smaller than the one in the normal phase when q < qa 
(see the right plot of figure 9). Take q = 1.150 < qa as an example, the grand potential 
presented in figure 9 shows clearly that the value of Q has a sudden jump from the normal 
phase to the condensed phase at u«a œ 1.8597, indicating a zeroth order phase transition. 

A phase diagram in terms of u and q for the case with m? = 0 is shown in figure 10. 
Indeed, for large enough q, the critical value of qu approaches to a constant, which is 
consistent with the probe limit shown in appendix D. This phase diagram is reminiscent 
of the case in the SU(2) holographic p-wave superconductor/insulator model presented in 
figure 12 in ref. [32]. Comparing our model with the special case m? = 0 to the SU(2) 
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Figure 8: The condensate of the vector operator in the case of m? = 0, q = 1.500 (left 
plot) and q = 1.200 (right plot). The left plot shows a second order phase transition, while 
the right plot shows a first order phase transition from the normal phase to the condensed 
phase. In the right figure, only the blue part is physical. 


Figure 9: The grand potentials with respect to the chemical potential in the case of m? = 0, 
q = 1.185 € (qa, %œ) (left plot) and q = 1.150 < qa (right plot). The grand potentials of 
the pure AdS soliton and hairy soliton are denoted by the green dashed line and blue solid 
line, respectively. In the case of q = 1.185, the discontinuity of the grand potential appears 
above the normal phase line. In the right plot for the case of q = 1.150, one discontinuous 
point at He = 1.8597 is under the normal phase line 2K?Q/V3 = —1. 
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Figure 10: The phase diagram for the case of m? = 0. The red, green and black lines 
stand for the second order, first order and zeroth order phase transitions, respectively.. 
Numerical results show uq ~ 2.2653 when 1/q approaches to zero. 


model, we find that the equations of motion in our setup look very similar with each 
other and the asymptotical expansions near the boundary are the same. Actually, we have 
carefully checked our numerical calculation in ref. [32] and found that similar discontinuous 
points in the grand potential also appear for the back reaction larger than a certain value. 
More precisely, the discontinuity first appears with grand potential larger than the normal 
phase, thus can not change the phase structure. However, for very large back reaction, 
i.e. larger value of a in ref. [32], the lower discontinuous point enters the region which 
has a lower grand potential than the one in the insulating phase. Thus we can also find a 
zeroth order transition from the normal phase to the condensed phase in the SU(2) p-wave 
model. If we add such a zeroth order transition region to the phase diagram in the SU(2) 
p-wave model, the figure 12 in ref. [32] is quite similar to figure 10 in the present paper. 
This result is consistent with our probe analysis in the presence of a magnetic field that 
the Einstein-Maxwell-complex vector field model is a generalization of the SU(2) p-wave 
model in the sense that the vector field has general m? [30]. 


3.3.3 Phase transition for m? < m?, 


When m? is less than m2,, there are two m?-dependent critical values of charge qe and qa 
(de < qa). The possible phase transitions and their orders are shown in table 3. For the 
above two cases with m? > m2, the condensed phase seems to survive even to the limit 
u — oo. ® In contrast, in the case with m? < m?2,, the hairy solution can not exist above a 
finite chemical potential. For the weak back reaction with q > qa, we first meet a second 


8Due to the precision of numerical control, it becomes much difficult to solve the system numerically 
for very large chemical potentials. Nevertheless, we have scanned a wide range of chemical potential up 
to 10° and found no evidence that the curves of grand potential and condensate would turn back to the 
small chemical potential region. 
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Phase transition and its order, m° < m2, 
a> fia, BC 


a 


<k 


Table 3: The phase transition and its order with respect to the charge and chemical 
potential in the case of m? < m2,. In the table, S=normal phase/insulator with pure 
AdS soliton solution, SC=condensed phase/superconductor with hairy soliton solution. 
0%, 1%, and 2 stand for the zeroth order, first order, and second order phase transition, 


respectively. 


order transition at He from the normal phase to the condensed phase, and then a zeroth 
order transition from the condensed phase back to the normal phase at uc2. In the range 
de < q < qa, the only difference is that the phase transition happening at ju; becomes first 
order. When we decrease the value of q less than qe, the condensed phase only appears at 
chemical potentials smaller than a critical one uc and has grand potential larger than the 
normal phase. The latter means the normal phase is thermodynamical preferred. We take 
m? = —3/4 as a typical case in the following. 

In this case, the two critical q are qa ~ 1.0801 and qe ~ 0.9951. When q > qa, there are 
two critical chemical potentials denoted by uea œ 0.8690 and ue ~ 1.8650, respectively. 
As the chemical potential is increased to ue, a second order phase transition from the 
normal phase to the condensed phased appears. If one continues scanning the shooting 
parameters, there is a critical value uez at which the curve of condensate turns back to 
the small chemical potential region. We can see from figure 11 that above ue the hairy 
solution does not exist, hence, a zeroth order phase transition from the condensed phase 
to the normal phase happens at pez. 

With the value of q is decreased into the region qe < q < qa, the first phase transition 
at Hea becomes a first order one (see figures 11 and 12). The zeroth order transition at pes 
still survives in this case. Our numerical calculation indicates that the difference between 
Ha and uez decreases as q decreases. Therefore, the region of chemical potential where the 
condensed phase is thermodynamical preferred is reduced by increasing back reaction. At 
the critical value q = qe, He is equal to e2, thus there is no thermodynamically favored 
condensed phase. 

For sufficiently large back reaction with q < qe, one can see from figure 13 that hairy 
solutions only appear below a critical chemical potential and their grand potentials are 
larger than those of the normal phase. Therefore in this case there is no physical phase 
transition from the normal phase to the condensed phase when q < qe. This phenomenon is 
an analog to the “retrograde condensation” appearing in the case of black hole background 
in ref. [28]. 

We construct the (u, q) phase diagram for the case with m? = —3/4 in figure 14. We can 
see from the phase diagram that when q > qe, there exist the insulator /superconductor/insulator 
phase transitions as one simply increases the chemical potential. For very large q, the back 
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Figure 11: The grand potentials with restive to the chemical potential in the case of 
m? = —3/4, q = 2.000 > qa (left plot) and q = 1.040 € (qe, qa)(right plot). The grand 
potentials of the pure AdS soliton and hairy soliton are denoted by the green dashed line 
and blue solid line, respectively. In the case of q = 2.000 > qa, there is a second order phase 
transition from the normal phase to the condensed phase at ueı œ 0.8690 and a zeroth 
order phase transition from the condensed phase to the normal phase at ue ~ 1.8650 . 
In the case of q = 1.040 € (qe, qa), there is a first order phase transition from the normal 
phase to the condensed phase at jie; ~ 1.6444 and a zeroth order phase transition from 
the condensed phase to the normal phase at ue2 ~ 1.7260 . 
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Figure 12: The condensation of the vector operator in the case of m? = —3/4, q = 2.000 
(left plot) and q = 1.040 (right plot). In each plot, the thermodynamically preferred hairy 
phase is drawn by the blue solid line. The left plot shows that the condensate arises from 
zero as u > 0.8690 and jumps to zero again when u > 1.8650, which represents a second 
order phase transition and a zeroth order phase transition happening at these two chemical 
potentials, respectively. In the right plot, the condensate has a jump at u ~ 1.6444 and 
u ~œ 1.7260, respectively, which corresponds to a first order phase transition and a zeroth 
order phase transition. 
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Figure 13: The grand potentials with respect to the chemical potential in the case of 
m? = —3/4, q = qe ~ 0.9951 (left plot) and q = 0.800 < qı (right plot). The grand 
potentials of the pure AdS soliton and hairy soliton are described by the green dashed line 
and blue solid line, respectively. In the case of q = qe, the two critical chemical potentials 
for the first order phase and the zeroth order phase transition coincides. In the case of 
q = 0.800, the grand potential of the condensed phase in the entire region is larger than 
the one of the normal phase. 


21 


Phase diagram (m?=-3/4) 
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Figure 14: The phase diagram for the case with m? = —3/4. The red, green and black 
lines stand for the second order, first order and zero order phase transitions, respectively. 
Numerical results show ueq ~ 1.7382 when 1/q approaches to zero. 


reaction of matter fields to the geometry is small, one should return to the results obtained 
from the probe limit. Indeed, we find that the critical chemical potential of the second 
order transition is very closed to a constant which is just the one by the probe calcula- 
tion in appendix D. For a very large q, the black curve corresponding to the zeroth order 
transition from the condensed phase to the normal phase still exists. But we can not see 
this critical point in the probe limit. Actually, our numerical results for q >> 1 suggest the 
following two scaling relations hold for two critical chemical potentials 


1 
Mer & =, fleg x gS, C(m?) = 0.4618m? + 0.0524m* + O(m®). (20) 
q 


with negative m?. The first relation for He is consistent with the probe approximation by 
using general Heun function analytically. The second relation for uez is found by numerical 
analysis. For more details please see appendix D. 

This scaling relation on [4.2 shows that the value of ue will be divergent for q —> oo 
once the value of m? is non-negative. That is to say, there does not exist similar turning 
point to the normal phase for m? > 0. This suggests that the value of m2, is indeed 
zero. Furthermore, for the case with m? < m2,, one can see that que x gS*+ is divergent 
as q — oo, which means that we can not see the turning point in the probe limit. The 
“retrograde condensation” can only appear for strong back reaction case. 


4 AdS Black Hole with Vector Hair 


4.1 Equations of motion and boundary conditions 


The hairy AdS charged black holes in the Einstein-Maxwell-complex vector field theory 
have been well studied in four dimensions in the probe limit [25] as well as in the case with 


22 


full back reaction [28]. Especially for the latter case, depending on two model parameters 
m? and q, the model manifests a rich phase structure. In this section we will construct 
asymptotically AdS charged black holes with non-trivial vector hair in (4+1) dimensions. 


Similar to ref. [28], we adopt the following metric ansatz 


d 2 

ds? = —a(r)e dt? + L 

a(r) 

The position of horizon is denoted as rp, at which a(r) = 0 and the conformal boundary is 

located at r + oo. For matter fields, we turn on the x component of p, and the temporal 
component of U(1) gauge field, i.e., 


+ r7(e(r)dz? + dy? + dz”). (21) 


Pdr” =p lr de, A, dz” = ġ(r)dt. (22) 


Without loss of generality, we can take py to be real. Then, the independent equations of 
motion are deduced as follows 
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where the prime denotes the derivative with respect to r. 

Those above equations of motion are very similar to those in the (3+1) dimensional 
spacetime [28]. Therefore, we use the same procedure in ref. [28] to solve the set of 
equations (23) by shooting method. The asymptotical expansion for metric fields and 
matter fields near the boundary r — oo is as follows 


P Pr-  Pre4 
=u- + er, Px = AS A Je , 
a4 ' C4 . Ki (24) 
asr’ +7) te cadh a Tor beoe i s 


where the dots stand for the higher order terms of 1/r. We choose the source free condition 
Px- = 0 as before. The coefficients u, p, and pr, are directly related to the chemical 
potential, charge density and x component of the vacuum expectation value of the vector 
operator Je in the dual system, respectively. 

We focus on black hole configurations that have a regular event horizon located at 
r, and require the regularity conditions at the horizon r = rp, which means that all 
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five functions {p,(r),(r), a(r), b(r), c(r)} will have finite values at r, and admit a se- 
ries expansion in terms of (r — ra). We have six independent parameters at the horizon 
{rn Pelra), (Tr), (rn), (72) }. Other coefficients can be expressed in terms of those pa- 
rameters. 

Two free parameters b(r;,) and c(r;,) can be fixed by AdS boundary conditions that 
b(r + oo) = 0 and c(r > co) = 1. Without loss of generality, the location of r, can 
be fixed to be one in our numerical calculation. We are then left with two independent 
parameters {pz(ra), ¢'(7,)}. By choosing @'(r;,) as the shooting parameter to match the 
source free condition, i.e., pz. = 0, we finally have a one-parameter family of solutions 
labeled by the value of pẹ at the horizon. After solving the set of equations, we can 
read off the condensate (J), chemical potential u and charge density p directly from the 
asymptotical expansion (24). 


4.2 Thermodynamics 


We will discuss how to extract thermodynamic quantities from our solutions. We will work 
in the grand canonical ensemble with chemical potential fixed. According to gauge/gravity 
duality, the Hawking temperature of the black hole is identified with the temperature of 
boundary thermal state, which is given by 


al (rp, )e7 brn) /? 


T= ; 
At 


(25) 
and the thermal entropy S is given by the Bekenstein-Hawking entropy of the black hole 


S=— A= EAS CTs | (26) 


4T ATV; 4 
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with A denoting the area of the horizon and V3 = f dxdydz. 

In order to find thermodynamically favored phase, we should calculate the grand po- 
tential Q in the grand canonical ensemble. Q of the boundary thermal state is identified 
with temperature T times the on-shell bulk action in Euclidean signature. The Euclidean 
action can be calculated by 


= Ik Spuctidean = | Pevar +12 + Lm) + f dry —h(2K — 6), (27) 


T—Co 


where h is the determinant of the induced metric on the boundary, and K is the trace of 
the extrinsic curvature. We have introduced the Gibbons-Hawking boundary term for a 
well-defined Dirichlet variational problem. Once again, the potential surface counterterm 
for the matter fields has been neglected here since this term does not make any contribution 
in our setup of the source free condition. 
By using the equations of motion (23), the on-shell action reduces to 
— 2K? Shel = BBV sac = 2r = Valoo, (28) 


Euclidean 
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with 8 = 1/T. Substituting the asymptotical expansion (24) into (28), we obtain 


y 
O = TSh = (U + ca — ba). (29) 


In the case with m? > —1, we find that c4 — b4 = 0 (see appendix C). In addition, for 


the normal phase (conductor) described by the RN-AdsS black hole solution in (5), one has 
a, = -r$ — aut and c4 = b4 = 0. 

In general, different solutions obtained by shooting method with ra, = 1 will have 
different chemical potentials. Since we work in the grand canonical ensemble, the chemical 
potential for each solution must be the same. Let us note that the system admits a useful 


scaling symmetry 


r—Ar, (t,2,y,z2) = A an) (¢,P2) >Ap) 2 > a, (30) 
with À a real positive constant. Under this symmetry, the revelent quantities transform as 
(T, p, Q) > A(T, 4,2), S78, pp, poy aA pay. (31) 


Therefore, we can set the chemical potentials of all black hole solutions equal by using 
above transformation rule directly. Without loss of generality, we shall fix the chemical 
potential to be u = 1. Finally let us mention that apart from the temperature T and 
thermal entropy S, this scaling transformation (31) is also true for the soliton case. 


4.3 Phase transition in black hole backgrounds 


Our numerical calculation shows that the phase structure of the model in the (4+1) dimen- 
sional black hole geometry looks precisely the same as the case in the (3+1) dimensional 
black hole spacetime [28]. Therefore in this section we will just give main results and some 
relevant details can be found in ref.|28]. The phase transition in the black hole background 
depends on the values of the mass square m? and charge q. There exists a particular value 
of m? denoted by m2. For m? is above or below m?2,, our system exhibits distinguished 
behaviors. Our numerical analysis suggests that m2, = 0, up to a numerical error, which 
is same as the value obtained in the (3+1) dimensional case [28]. For m? > m3, the 
condensed phase exists even for sufficiently low temperatures, while for m? < m?,, the 
condensed phase will be absent when the temperature is lowered to a critical value. All 
possible phase transitions are summarized in table 4 and table 5. 

In the case of m? > m3, there is a mass dependent critical charge qe. When the charge 
of the vector field is larger than qe, the condensed phase will become thermodynamically 
preferred when the temperature is less than a critical temperature at which a second 
order phase transition occurs. When we decrease q to a certain value smaller than qe, 
the transition from the normal phase/conductor to the condensed phase/superconductor 
becomes a first order one. For large m? > m?,, the condensed phase seems to survive 
even if the temperature goes to zero, T — 0. In contrast, in the case with m? < m3, 
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Phase transition and its order, m2, < m? 
T>T., BHT =T?" T<T., BC 


Cue <a [P>T, BA[T=T.1] T<Te, BO 


Table 4: The phase transition and its order with respect to the charge and temperature 
in the case of m2, < m°. In the table, BH=normal phase/conductor with the RN-AdS 
black hole solution, BC=condensed phase/superconductor with hairy black hole solution. 
1% and 2” stand for the first order and second order phase transitions, respectively. 


Phase transition and its order, m° < m3 
T > Ta, BH | T= To,2" | To <T < Too, BC T < To BH 


ESR EES PaA 


q < dg 


Table 5: The phase transition and its order with respect to the charge and temperature 
in the case of m2, > m?. In the table, BH=normal phase/conductor with the RN-AdS 
black hole solution, BC=condensed phase/superconductor with hairy black hole solution. 
0%, 1% and 2” stand for the zeroth order, first order and second order phase transitions, 
respectively. 


the condensed phase cannot exist when the temperature is lower than a certain value. To 
determine the precise value for m2, , we need to solve the coupled equations of motion (23) 
at very low temperatures to see whether the condensate would turn back to the higher 
temperature region. The precision of numerical calculation prevents the investigation on 
the limit of T — 0. Nevertheless, our numerical calculation suggests that m2, is very 
closed to zero. We will consider one concrete example for both cases with m? > m2, and 
m? < m?2,, respectively. 


4.3.1 Phase transition for m? > m3 


In the case with m? > m?,, as we lower the temperature, we can always find a phase 
transition from the normal phase to a condensed phase where a charged black hole develops 
non-trivial vector hair and is thermodynamically favored. In the dual field theory side, 
it means that a vector operator acquires a vacuum expectation value (Je) # 0 which 
breaks the U(1) symmetry as well as the rotation symmetry in x — y plane spontaneously. 
Furthermore, the order of the transition depends on the strength of back reaction controlled 
by q. For small back reaction, this transition is second order, while for large back reaction, 
it becomes a first order one. The critical value of q denoted as qe depends on m? we choose. 
We here consider m? = 5/4 as a concrete example, where qe ~ 2.0972. One can find similar 
results for other values of m? > Gee 

Let us first consider, for example, the case with q = 2.500 > qe. The condensate appears 
at Te ~ 0.0698 and arises gradually from zero for lower temperatures. Comparing grand 
potential between the hairy solution and the RN-AdS solution, one can see from the left 
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Figure 15: The grand potential (left plot) and condensate (right plot) as a function of 
temperature in the case of m? = 5, and q = 2.500. In the left plot, the grand potential 
of the RN-AdS black hole solution and hairy black hole solution are described by the 
green dashed line and blue solid line, respectively, and the phase transition happens at 


Te ~ 0.0698. 


plot of figure 15 that the condensed phase is thermodynamically favored. It is a second 
order phase transition with the critical behavior (J;) ~ (1 — T/T,)"/2 near T,. The critical 
exponent for the condensate is 1/2, the same as the one in mean field theory. 

A qualitative change happens as q is decreased to below qe. The grand potential and 
condensate as a function of temperature for q = 1.800 < q are shown in figure 16. We can 
see that a swallowtail appears in the grand potential and in this case, the condensate has 
a jump from zero to a finite value at the critical temperature, which is a typical feature of 
a first order phase transition. This tells us that when q < qe, the phase transition becomes 
first order from the normal phase to the condensed phase as the temperature is lowered. 

A phase diagram in terms of temperature and charge is shown in figure 17 in the case 
with m? = 5/4, where the red and blue curves stand for the second order and first order 
phase transitions, respectively. 


4.3.2 Phase transition for m? < må, 


The phase structure for m? < mê is the same as the one in the (3+1) dimensional black 
hole background [28]. The parameter space of q is divided into three regions by two special 
values of q, denoted as qa and qg with qa > qg. The phase behavior changes qualitatively 
in each region. To compare our results to the case with one dimension less, we take the 
same terminology, i.e., the critical temperatures for second order transition, first order 
transition and zeroth order transition are denoted by Tes, Te and To, respectively. In this 
subsection, we consider m? = —3/4 as a concrete example. In this case, qa ~ 1.766 and 
gg ~ 1.700. 
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Figure 16: The grand potential (left plot) and condensate (right plot) as a function of 
temperature in the case of m? = 5, and q = 1.800. In the left plot, the grand potential of 
the RN-AdS black hole solution and hairy black hole solution are described by the green 
dashed line and blue solid line, respectively. The phase transition happens at T ~ 0.0145. 
The physical branch for the condensate in the right figure is denoted by the blue solid line. 
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Figure 17: The phase diagram for the case with m? = 5/4. The red line and blue line 
stand for second order and first order phase transitions, respectively. When 1/q — 0, we 
have T./q œ 0.0613. 
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Figure 18: The grand potential (left plot) and condensate (right plot) as a function of 
temperature in the case of m? = —3 and q = 2.000. In the left plot, the grand potential of 
the pure RN-AdS solution and hairy solution are described by the green dashed line and 
blue solid line, respectively. The condensate begins at Ty. = 0.1592 with a second order 
phase transition and ends at Teo = 0.1362 with a zeroth order phase transition. In the 
right figure, only the blue solid line stands for the physical branch of the condensate. 


For the weak back reaction q > qa, we show the condensate and grand potential versus 
temperature in figure 18. We immediately see that the condensate arises below Ty, and 
then turns back to the high temperature region at Tọ (see the dashed orange dashed line, 
which corresponds to the upper branch of the grand potential in the condensed phase). 
When T < To, there does not exist hairy black hole solution and thus in the region with 
T < T.9 we have only the RN-AdS black hole solution. By comparing the grand potential 
for each solution, we find that the thermodynamically favored region of the condensed phase 
only locates in the region Tig < T < Te. At Tyo, it is a second order phase transition, 
while there is a zeroth order transition at T.9, since at the critical temperature the grand 
potential has a sudden jump from the condensed phase to the normal phase. 

As we decrease the value of q to the region q8 < q < qa, there is a little difference 
compared to the previous case. As we lower the temperature, we will first see a first order 
transition from the normal phase to the condensed phase at 7’.,, and then a zeroth order 
transition back to the normal phase at To. A typical case with q = 1.719 is shown in 
figure 19. 

In the above two cases, as we lower the temperature, there exists a condensed phase 
which is thermodynamically favored, although only in a narrow temperature range. The 
story has a dramatic change as q is decreased past qg. In this case, the hairy black hole 
solution only appears in the high temperature region. In figure 20 we present the grand 
potential as a function of temperature in the cases of q = 1.700 (left plot) and q = 1.600 for 
m = —3. One can clearly see that the grand potential of the hairy black hole solutions is 
always larger than the one for the RN-AdsS black hole at the same temperature. Thus these 
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Figure 19: The grand potential (left plot) and condensate (right plot) as a function of 
temperature in the case of m? = —3 and q = 1.719. In the left plot, the grand potential 
of the pure RN-AdS solution and hairy solution are described by the green dashed line 
and blue solid line, respectively. The condensate begins at Tea = 0.1194 with a first order 
phase transition and ends at Teo = 0.11856 with a zeroth order phase transition. In the 
right plot, the solid blue line is the physical branch of the condensate, while two dashed 
orange parts are un-physical. 


hairy black hole solutions are thermodynamically disfavored. This is just the so called 
“retrograde condensation” reported in ref. [28]. Let us notice that similar “retrograde 
condensation” also happens in the soliton case in the sense that the hairy soliton solution 
only appears below a critical chemical potential and is thermodynamically disfavored, as 
shown in figure 13. 

A phase diagram in terms of temperature and charge is shown in figure 21 in the case 
with m? = —3/4. There the red, blue and black curves represent second order, first order 
and zeroth order phase transitions, respectively. 


5 The Complete Phase Diagram 


The above studies tell us that, depending on two model parameters m? and q, the Einstein- 
Maxwell-complex vector field theory exhibits rich phase structures. From the perspective 
of dual field theory, the system admits four phases including an insulating phase described 
by the pure AdS soliton, a soliton superconducting phase described by the hairy soliton, a 
conducting phase described by the RN-AdS black hole, and a black hole superconducting 
phase described by the hairy black hole. For given m? and q, the complete phase diagram 
can be constructed in terms of temperature and chemical potential. At each point in T-u 
plane, one should find the phase which has the lowest grand potential. Since the dual 
system in spatial directions is homogeneous and infinite, the total grand potential Q is 
divergent. So we consider the grand potential density defined by w = 2/V3 for both the 
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Figure 20: The grand potentials as a function of temperature in the case of m? = —3, 


q = 1.700(left) and q = 1.600(right). In both figures, the grand potential of the RN-AdS 
black hole solution and hairy black hole solution are denoted by the green dashed line and 
blue solid line, respectively. 
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Figure 21: The phase diagram for the case with m? = —3/4. The red and blue lines stand 
for second order and first order phase transitions, respectively, while the black one for the 
zeroth phase transition. When 1/q — 0, we have Ty.2/q ~ 0.1114 and T.9/q > 0. 
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soliton and black hole cases. Remember that we fix the size of the compact directions 7 
and z to be T= 7. 

As a warmup, let us first show some properties of the grand potential density w(T, u). 
We note that one may associate an arbitrary temperature to the Euclidean soliton for a 
given chemical potential. So the grand potential density of the soliton case is independent 
of temperature, i.e., Wsoliton(T, U) = Wsolitonl u). Therefore, the phase boundary between 
insulator and soliton superconductor should be a line parallel to the T axis. According 
to the scaling law (31), the grand potential density of black hole for different chemical 
potential is given by wga (T, u) = wgy(T, u = 1)u*. Taking advantage of above properties, 
one can easily find a first order phase transition between the AdS soliton and RN-AdS 
black hole, °? and the phase boundary between the RN-AdS black hole and the AdS soliton 
at zero temperature is at a chemical potential u = 121/4 ~ 1.8612, while the boundary 
at zero chemical potential locates at T = 1/7. In addition, note that since the critical 
temperature between the RN-AdS black hole and the hairy black hole is proportional to 
Lt, this phase boundary between them is a straight line passing through the original point 
(T, u) = (0, 0). 

Now we begin to construct the T-u phase diagram. To find the thermodynamically 
favored phase among four phases, we adopt a method similar to the one in ref. [24]. In the 
first step, we scale all solutions to the ones with the same value of chemical potential by 
using above scaling rules. Then we scan temperature to find the physical phase boundary 
at this u by comparing grand potential densities for all solutions. Repeating this procedure 
with different values of u, we can assemble the full T-u phase diagram. 

Remember that for different choice of model parameters m? and q, the phase behaviors 
can qualitatively change. More precisely, for regions m? > m?2,, 0 < m? < m2, and m? < 0, 
our system exhibits distinguished thermodynamics. It is a very hard work to construct the 
T — u phase diagram for all values of m?. Instead, we consider some typical examples for 
each parameter range. We hope that the phase diagram will not have a dramatic change 
for other choice of m? in each region. We draw some typical examples in figure 22. 

In the region with m? > m?2,, let us consider m? = 5/4 as an example. The phase 
structure depends on the strength of back reaction. For weak back reaction with q > q2 ~ 
1.5345, the phase diagram looks like the top-left plot in figure 22. In this range of q, as we 
lower q, the phase boundary between the insulator and soliton superconductor goes toward 
larger u, and the slope of the line separating black hole superconductor from conductor is 
reduced. This is consistent with the fact that the back reaction makes the phase transition 
more difficult. Furthermore, the boundary between the two kinds of superconducting 
phases moves to a lower temperature for smaller q. As q is lowered past q2, there will be 
two discontinuous points in grand potential indicating zeroth order phase transition (see 
figure 3). Between two discontinuous points, there is no soliton superconducting phase, 
thus it changes the phase diagram a lot. As shown in the top-right plot of figure 22, the 
soliton superconducting phase is broken into two parts and there is a particular range of 


°This first order phase transition is the planar analogous of the Hawking-Page transition for black holes 
in global AdS coordinates [36]. 


32 


Phase diagram(m7=5/4, q=2) 


Phase diagram(m*=5/4, q=1.534) 


Phase diagram(m*=-3/4, q=2) 


re g? Phase diagram(m7=5/4, q=1.5) 05 
0.4} BH 
6- 
0.3! 
T| E 
0.2 
27 0.1} sc 
% % 05. i 15 2 25 3 
u 
Phase diagram(m*=-3/4, q=0.8) 
0.35} 
BH 
0.3 BH 
0.25+ 
T 0.2 
0.15+ 
0.1 
0.05+ s 
(0) L L J L ri 1 
0 0.5 1 1.5 2 2.5 2 2.5 


Figure 22: The complete phase diagram for the model with different mass square and 
charge of the vector field. In the figure, S=pure AdS soliton, BH=RN-AdS black hole, 
SC=hairy soliton, BC=hairy black hole. 
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chemical potential where the thermodynamically favored phase is replaced by black hole 
superconductor for lower T and conductor for higher T. When the value of q is decreased 
to the case that the transition from the pure AdS soliton to the hairy soliton becomes 
zeroth order (see figure 4), for some suitable values of charge in this region, the insulating 
phase and soliton superconducting phase separate from each other, between them there 
are phases coming from the black hole background. A typical example for this is shown in 
the middle-left plot in figure 22. 

In the region of m2, < m? < m?,, we take m? = 0 as an example. We find that in this 
case the T-u phase diagram is qualitatively similar to the one numerically constructed for 
the holographic s-wave model [24] as well as the one schematically plotted in ref. [33] for the 
SU(2) p-wave model. The four phases typically meet in two triple points. As we decrease 
q, the two triple points gradually approach to each other, merge into one quadruple point 
where four phases can coexist and then pass through each other. There exists a phase 
boundary in which the black hole superconductor becomes an insulator via a first order 
phase transition. 

In contrast to above two cases, when m? is negative, the hairy solution can exist only 
in a certain chemical potential region for soliton background and temperature region for 
black hole background. As one can see from figure 13 and figure 20 that the hairy solution 
becomes thermodynamically irrelevant for sufficiently strong back reaction, even more in 
the near zero temperature region. Due to this, the T-u phase diagram has a dramatic 
change compared to the cases with m? > 0. Let us consider the case m? = —3/4 as an 
example. Indeed, the structure of the phase diagram depends on the value of q we choose. 
A typical diagram for small back reaction is drawn in the middle-right plot of figure 22. 
One can see that the lower-right region of the phase diagram is replaced by the conducting 
phase. As we lower g, the region in which the two superconducting phases dominate is 
deduced. In the case with q = 1.150 shown in the bottom-left plot of figure 22, the soliton 
superconducting phase shrinks to a smaller region, while the black hole superconducting 
phase disappears. For sufficiently small q, both the black hole superconducting phase and 
the soliton superconducting phase become thermodynamically disfavored compared to the 
uncondensed phases. In this case, we are left with only the conducting phase and insulating 
phase at hand. There is a first order phase transition between them shown in the bottom- 
right plot in figure 22, which is the well-known confining/deconfining phase transition in 
gauge field theory [23]. 


6 Summary and Discussions 


In this paper we give a detailed study on the holographic p-wave superconductor model 
based on the Einstein-Maxwell-complex vector field theory with a negative cosmological 
constant. We work in five dimensional asymptotically AdS backgrounds with full back 
reaction of matter field on the background geometries. In this model, we find four kinds 
of bulk solutions given by the pure AdS soliton, RN-AdS black hole and their vector 
hairy counterparts. According to the AdS/CFT dictionary, the hairy solution is dual to 
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a system with a non-zero vacuum expectation value of the charged vector operator which 
breaks the U(1) symmetry and the spatial rotation symmetry spontaneously. The above 
four solutions in the bulk correspond to an insulating phase, a conducting phase, a soliton 
superconducting phase and a black hole superconducting phase, respectively. 

In this model there exist two model parameters, i.e., the charge q and the mass square 
m° of the charged vector filed p,. The phase structure of the model heavily depends on 
the two parameters. In the black hole background case, the phase behavior is qualitatively 
same as the one in the (3+1) dimensional case studied in [28]. One can see second order 
transition, first order transition, zeroth order transition as well as the “retrograde conden- 
sation”. The latter two cases can only appear for the case m? < m? with m? = 0 suggested 
by numerical analysis. 

For the soliton case, our system exhibits distinguished behavior depending on concrete 
value of m?. Qualitatively, the parameter space for m? is divided into three regions with 
m? > m2, = 0.218 + 0.001, 0 < m? < m?, and m? < 0, respectively. In each region, the 
phase behavior depends on the strength of back reaction g. There also exist second order, 
first order and zero order phase transitions as well as the “retrograde condensation”. An 
interesting observation is that when m? < 0, the hairy black hole solution does not exist 
as the temperature is lower than a certain value (see figure 18), which leads to the zeroth 
phase transition at the critical point. We draw the phase diagrams (see figures 6, 10 and 
14) in each region of mass parameter space. Compared to the black hole background case, 
an additional complication appears in the soliton background case. 

Recall that our one parameter family of hairy solutions are labeled by the value of pz 
at the tip of the soliton. For m? > 0, as we lower q to a critical value, there exists a region 
of chemical potential, in which the hairy soliton solution does not exist. As a result, the 
curve of the condensate (or grand potential) parameterized by p,(r9) breaks down, which 
results in two discontinuous points (see figures 3, 4, 5, and 9). For the existence of this 
discontinuous point and the hairy soliton solution at this point, we give a detail analysis 
in appendix A and appendix B. Some physical meanings of this discontinuous point both 
at the gravity side and the dual field theory should be further studied. For the case 
with m? < 0, the hairy soliton solution terminates at a finite chemical potential for weak 
back reaction, while for sufficiently large back reaction, it becomes thermodynamically 
subdominant compared to the pure AdS soliton. See figures 11 and 13 for details. 

With all four phases at hand, we construct some complete phase diagrams in terms 
of chemical potential and temperature. As one can see, there are many types of phase 
transitions in both soliton and black hole backgrounds, which depend on the values of 
parameters m? and q. The T- phase diagrams are expected to be much more complicated 
than the ones for the holographic s-wave model and SU(2) p-wave model [24, 33]. Some 
typical examples are shown in figure 22. It has been suggested that the Einstein-Maxwell- 
complex vector field model is a generalization of the SU(2) p-wave model in the sense that 
the vector field has a general mass and gyromagnetic ratio y. Comparing our complex 
vector field model to the SU(2) model with a constant non-Abelian magnetic field, we find 
that the SU(2) p-wave model can be recovered by the restriction m? = 0 and y = 1 at 
least in our setup in ref. [30]. Similar statement holds for the black hole case. Indeed, for 
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the special case m? = 0, our model exhibits very similar phase structures as the SU(2) 


model [32, 35]. In particular, the complete phase diagrams in our model with m? = 0 are 
qualitatively same as those drawn schematically in ref. [33]. We have freedom to choose 
other values of m? in our model, thus it can be used to describe much rich phenomena 
in the dual strong coupled systems. As one can see, depending on the values of m? and 
q, we have indeed many other types of phase diagrams which do not appear in the SU(2) 
p-wave model [33] as well as in holographic s-wave model [24]. In addition, let us mention 
that in the soliton background case, we have checked that there also exists a zeroth order 
phase transition in the SU(2) p-wave model as the back reaction is strong enough, which 
is missed in the previous studies. 

Finally some remarks are in order. (1) There exist two critical masses of the vector 
field in our model, m?, = 0.218 + 0.001 and m2, = 0, in the soliton background case, 
a critical mass, m?, = 0, in the black hole background case. Note that the latter also 
appears in the (3+1) dimensional black hole case [28], there some possible implications 
of the critical mass have been discussed and therefore we will not repeat here. While the 
critical mass m2, = 0.218 + 0.001 looks spacetime dimension dependent, its appearance is 
interesting and is worthy to further investigate. Unfortunately at the moment we have no 
any physical interpretation for this. (2) In both the soliton and black hole backgrounds, our 
model exhibits zeroth order phase transition, which has a discontinuity of grand potential 
(free energy) at the critical point. In ref. [38] the author argues that a zeroth order phase 
transition could appear in the theory of superfluidity and superconductivity and presents 
an exactly solvable model for such a phase transition. The authors of ref. [39] show that a 
zeroth order phase transition exists in an exactly solvable pairing model for superconductors 
with p, + ipy-wave symmetry. (3) In the soliton background case, our model shows the 
insulator /superconductor/insulator/superconductor phase transitions when the chemical 
potential continuously increases. Indeed such a series of phase transitions is expected to 
appear from theoretical perspective in some superconducting materials, for example, see 
the phase diagram in the figure 10 of ref. [40]. (4) We can see from the complete phase 
diagrams in terms of temperature and chemical potential that in some cases, more than one 
superconducting phase appear in a phase diagram in our model. The phase diagram for 
some realistic superconducting materials is usually complicated, and indeed, more than one 
superconducting phases can occur, for example, see ref. [41, 42, 43, 44, 45, 46]. Definitely, it 
is of great interest to see whether our model is relevant to these superconducting materials. 
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A The regular boundary conditions for the equations 
of motion 


In this appendix, some details for solving the coupled equations of motion (11) will be 
given. In particular, we will show why in the soliton background case, the hairy soliton 
solution does not exist in some region of chemical potential, shown in subsection 3.3.1. 
The black hole case will be compared in the end of this appendix. 

In order to find the solutions for all the six functions F = {g(r), f(r), h(r), x(7), pe(r), (r) $, 
one must impose suitable boundary conditions at boundary r —> oo and the tip r = rp. The 
asymptotic behaviors for these functions near the infinite boundary can be easily obtained 
from the equations (11), which are shown in equations (12). The condition of p,_ = 0 
leads the problem to be a boundary value problem. In this paper, the method to solve 
them is shooting method, i.e, by finding some suitable initial values at the tip for each 
function in F in order to give result of p,_ = 0 near the boundary when we integrate the 
equations (11) from the tip to the boundary. To solve the equations (11), one needs ten 
initial conditions, i.e, {g(ro), f (ro), h(ro), x(ro), Pz(ro), (70); f (ro), R (ro), ph(ro), O' (70) }- 

By the scaling symmetries in (14)-(17), one can first set {ro = 1, f(ro) = 1, h(ro) = 
1,x(ro) = 0} . The first order derivatives of f(r), h(r), px(r) and g(r) at the tip can be 
obtained by imposing the requirement that each term in equations (11) behaves regular at 
the tip. By some calculations, we have 


Fro) = $'(r0) bo = 2" popolaro), (32) 
h'(ro) = —2p'.oP0 = 2p29(m” — 4°¢5)/9'(r0), (33) 
g'(ro) = 4 —  02065/3, (34) 


where pro = Px(ro) and œo = ġ(ro), which are two free parameters at hand. We can choose 
fixing pro using (ro) as the shooting parameter or vice versa to match the source free 
condition, i.e., Px_ = 0. After solving the coupled differential equations, we should use 
the scaling symmetries (14)-(17) again to satisfy the asymptotic conditions f(oo) = 1, 
h(co) = 1 and y(co) = 0 and the fixed compacted length IT = z of the coordinate 7. Then 
we can obtain the condensate (J), chemical potential u, charge density p, and grand 
potential by reading off the corresponding coefficients in (12), respectively. 

The equation (34) should be considered carefully, as it gives a restriction on the solution 
for equations (11). By definition, the tip ro is the maximal zero point of g(r). Therefore, 
g'(ro) must be non negative. For the non-extremal case, we should have a positive value 
of g'(ro), which gives a restriction for the initial values of pvo and ġo, 


4 — q?p262/3 > 0 (35) 


It is precisely this restriction so that we cannot get the hairy soliton solution within some 
region of chemical potential in the condensed phase, which is shown in figures 3, 4 and 9. In 
order to show this discontinuous phenomenon more visualizable, in figure 23 we show two 
kinds of orbits for the initial values p,9 and ¢o for the case with m? = 5/4 and q = 1.900 
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Figure 23: The orbit of the initial values pro and ġo for the shooting method in the case 
of m? = 5/4 with q = 1.900 (left plot) and q = 1.520 (right plot). The blue solid lines are 
the orbits for the initial values pzo and ġo, and the red dashed lines stand for the curves 
g'(ro) = 0. The regions below the red dashed lines give g'(ro) > 0. In the left plot, the blue 
line is continuous and has no any intersection with the curve g'(ro) = 0. In the right plot, 
the blue lines have two intersection points at @9 ~ 1.5378 and o ~ 1.3025, respectively. 


and 1.520, respectively. In the figure, The blue solid line is the orbit for the initial values 
Pro and ġo, while the red dashed lines stands for the curve g'(r9) = 0. The region below 
the red line has g‘(ro) > 0. 

From the figure 23, one can see that the orbit of ġo and pro behaves well and has no any 
intersection with the dashed red line when q = 1.900. In this case, the solution orbit leads 
to the continuous grand potential and condensate in the condensed phase, as shown, for 
example, in figure 1. However, when q = 1.520 as shown in the right plot of figure 23, the 
solution orbit is broken into two parts by the curve g'(ro) = 0. It leads to the discontinuity 
of the grand potential in the condensed phase, just as shown in figures 3, 4 and 9. Between 
these two intersection points, there do not exist any hairy soliton solutions and thus in 
that region we have to replace hairy soliton solutions by the pure AdS soliton because the 
latter can associated with any chemical potential. 

For the black hole case, we can also construct similar boundary conditions at the black 
hole horizon and AdS boundary. In this case, we can use the values of p,(rp_,) and ġ'(rp) 
as the shooting parameters. There is also a restriction at the horizon as 


a(t) = 4— LAr) >0. (36) 


One can see this restriction only depends on ¢’(r;,), which is very different from the one 
(35) in the soliton case. The numerical results in the region where we have scanned show 
this restriction does not cause any discontinuity in the condensed phase in the black hole 
case. 


38 


25 


20 Ing’ =- In(r-1,) 


ae 


decreas g’ C) 


—40 35 =30 -25 -15 -10 -5 


-20 
In@—-1p) 
Figure 24: The asymptotic behavior of x’. From the figure, one can see that when g'(ro) > 
0, there is an asymptotic behavior x’ = ae Here the values of g’(r9) are chosen from 
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B The extremal hairy soliton solution 


In the main text we have assumed g'(ro) > 0 in seeking for the solutions with the regular 
condition at the tip. In this appendix we will give the extremal solution with g'(ro) = 0. In 
order to find the extremal solution, it is helpful to observe the behaviors of these functions 
in F when g'(ro) approaches to zero. 

In figure 24 we plot the behavior of lIn y’ near the tip when g'(ro) tends to zero. It 
shows that when g'(ro) goes to zero, the behavior of x’ behaves as x’ = T Taking this 
behavior into account, we propose the following expansion near the tip 


gr) =a) =n) bese Ar = ro) Ee 
Plr) = po + bio(r — ro) 2 +++, Pol?) = pro + Prro(r = ro) +, (37) 
f =1+ fiolr—ro)*? +, h(r) =1+ Aio(r = ro 7 +- 


Here the scaling symmetry has been used to scale f(r) and h(r) to be unity at the tip. 


Substituting these assumptions into the equations (11), we find only 7 = 4 can give a 


2 
self-consistent expansion with the coefficients 
4pxo(m* — 0) 


89202 
10 = fio/d0 = ete Px10 = —hio/2p.0 a an 
i 93 (38) 


4 
O00 — 12 =0, gj = re / 489? — 24m? — 4m?q?¢2 — 720? + 3q465 + mt. 


In this case with j = 1/2, a remarkable fact is that although g'(ro) = 0, the solution with 
the expansion (37) still has a finite compacted length I = 87/(g;ro) for the coordinate n. 
In addition, the solution looks singular at the tip because y(r) is divergent there. In fact, 
this is a coordinate singularity and the solution is regular at the tip. To see this, one can 
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1/4 near the tip r = ro, the metric then has the following form 


introduce € = (r — ro) 
16 
ds? = ae + r2(—dé? + dz? + dy? + g;€7dn’). (39) 
J 

One can use the expansion (37) with the coefficients (38) to solve the equations (11) 
with the source free condition by scanning the shooting parameter ġo. 1? In the case of 
2 = 5/4 and q = 1.520, we find the first two values of o are 1.302900 and 1.537516, 
which are very closed to the corresponding values in figure 23, the latter are obtained by 

setting g'(ro) = 4 x 1074. 


C The stress energy tensor 


In this appendix, we study the stress energy tensor of the dual field theory on the boundary. 
From the AdS/CFT correspondence, the stress energy tensor can be obtained from the 
Brown- York tensor of the bulk spacetime. In the asymptotic AdS; spacetime, the Brown- 
York tensor on the boundary OM can be computed as [47] 


1 = 3 1 
[LV [LV he” pV pw : 4 
DRR =z(K -K +—h + 5G ) (40) 


Here ht” is the boundary metric, GH” is the Einstein tensor for the boundary metric h””. 
K” is extrinsic curvature, defined by the outward pointing normal vector to the boundary 
n as 


K” = Lva + VÝR”). (41) 
In our case, the boundary spacetime is flat, so the Einstein tensor for the boundary metric 
vanishes. In what follows we discuss the stress energy tensor for the soliton and black hole 
backgrounds separately. 
The soliton case 
By the metric ansatz (9) and the asymptotic expansion for the metric in (12), the 
leading terms of the Brown-York tensor turn out to be 


2K? T” = 7 §(Ahy + g4 — ya : 

2K T” = nn — 4 

2K? T”! = rÊ (4%4 — EA h4 
2K? T™ = oni — 4h4 — 


Fi 


) 
i: (42) 
). 


In the following, we omit the factor 2x? for convenience. The non vanishing components 


10To avoid the divergence of y’, f’, h’, ¢', p' at the tip in the numerical calculations, it is better to 
change the coordinate from r to € and make a suitable variable substitution on y such as y — eX/4. 
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of the stress energy tensor t” for the dual field theory on the boundary are 


t” = ha + g4 — 4X4, 
t” = Aya — ga — 4 fa, 


t = 4X4 — ga — 4 fa — 4ha, A 
t™ = 394 — 4h4 — 4 fa. 
The trace of the boundary stress energy tensor is 
Trt) = 12(x4 — fa — ha). (44) 
Note that the equation of motion for y(r) reads 
yf 2g! Woe” alto? | 8 8 _ (45) 


T g h 3rh 3rf a fa r 


Substituting the expansion (12) into the above equation, in the limit r => oo, we have 


2,2 
A(fa + hy — x4) = EE Pe 4 (4), 6 = VTE we (46) 
It is clear from the above equation that the dual field theory is traceless as 6 > 0. When 
6 = 0, the trace looks non-vanishing with Tr(t"”) = 22. But this is not true. The reason 
is as follows. Note that the case 6 = 0 is equivalent to the lower bound with m? = —1. 
As mentioned above, in this case, the expansion in (12) is no longer valid, instead some 
logarithmic terms should appear in the asymptotic expansion. In fact, it can be shown 
that with these additional terms the dual stress energy tensor is still traceless. 

The black hole case 

In this case, by the metric ansatz (21) and the asymptotic form of the metric in (24), 
we have the nontrivial terms of the Brown-York tensor near the boundary 


T* =r *(4e, — 3a4 
T*® = r (4b, — a4 
TY =r [4 (b4 — c4) — a4 
TA = r84 (ba — c4) — a4 


j, 
(47) 
i, 


where we have omitted the factor 2x? for convenience. The stress energy tensor t!” for the 
dual field theory on the boundary reads 


t= 4c, — 3a4, 
t” = 4b, — aa, 
t = [4(b4 — ca) — aul], 
t = [4(b4 — c4) — aa), 


(48) 


Al 


with its trace 


In this case, let us consider the equation of motion for the metric function b(r) 


Qa’ d 2p? re? = De’g? ppe 8r 4 

Ee ee ee (50) 
a c “are 3a 3ra2c a r 
Substituting the expansion (24) into the above equation, and taking the limit r — oo, we 
have the following asymptotic form 


2 +0) Pox 4 iby =V F m. (51) 


4(c4 = b4) =l 375 © 


As the soliton case, one can see that the stress energy tensor is also traceless in the black 
hole background case. 

Finally let us mention an interesting fact on the stress energy tensor of the dual field the- 
ory. Naively thinking, one may expect that the stress energy tensor should be anisotropic 
with t?” Æ t” in both hairy soliton and black hole cases, for which the metric functions h(r) 
and c(r) are nontrivial and from the point of view of dual field theory, the x-component of 
the vacuum expectation value of the vector operator is non zero, while the y-component 
vanishes. But this anisotropy does not appear in the stress energy tensor. Our numerical 
results suggest that h4 and c4 are zero, !! which leads to t®* = ¢¥ in both cases. 

A similar phenomenon has been reported in the paper [48], the authors have given an 
expression for the variation of the free energy of periodic AdS black holes with respect to 
the periods of transverse coordinates, which gives an explanation for this phenomenon in 
the black hole and soliton backgrounds. 


D The scaling relations between the critical chemical 
potential and charge 


In this appendix, we will give some discussions on the scaling relation (20), which holds for 
the large q case. We will focus on the soliton case. The relation between ue and q can be 
obtained by semi-analytical method through the general Heun function, while the relation 
between ue and q can only be obtained by numerical method. We first study the former 
case. 

In the case with large q, the probe limit is a good description near the phase transition 
point. In that case, the chemical potential takes the critical value, while the vector field is 
very small, whose influence on the background and electromagnetic field can be neglected. 


Tn the soliton case with m? = —3/4,q = 2.000 , the numerical results show that f4, g4 and y4 are of 
order one, but h4 is zero up to a numerical error 10712, while in the black hole case with the same model 
parameters, the numerical results show that a4 is of order one, but c4 and b4 vanish up to a numerical 
error 10712., 
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This means that we can just treat the vector field as a perturbation and take the background 
solution as 


1 
gr) =1— 4, fr) =Alr) = 1, x(r) = 9, g(r) = pa. (52) 
Then the equation for the vector reads 
/ 2,2 2 
n g Df q Her m 
2 ae ~~ =p, = 0, 53 
Pale ahah agar rag? (53) 
By introducing 
1 
r=1, pelt) = F(z), 6= VIF me, (54) 
2 


the equation (53) can be rewritten as 


z4(5 +28) — (26 + 1) 


27(6 +1)(6 +3) — uA 
eee Sey t eap 
z(2z4 — 1) 


p" 


F'(z)+ (z)=0. (55) 


This equation has a general solution as 


1 6+3 6+1 
F(2) =CyHeunG(—1, -79" 2), >, =, 5 + 1,1; =) 


4 
1 =g i= 
POr HeunG(—1, a a 1 - 6,1; -2”). 


(56) 


Here C and C are two constants and HeunG stands for the general Heun function with 
the feature HeunG(, , , , , ;0) =1. To satisfy the source free condition at the boundary, 
it is easy to see that one has to take Cj = 0. On the other hand, in order to have a 
regularity for the vector field pz at the tip, there exists a constraint condition at z = 1 as 


[24(5 + 26) — (25 + 1)] F"(z) — [2°(6 + 1)(6 + 8) — ual F(2)|2-1 = 0. (57) 


Substituting the solution (56) with C2 = 0 into (57), we can obtain a nonlinear equation 
for y2,q°. The equation can be easily solved by numerical method. The solution of course 
depends on 6 or m?. As a result we can write the solution as 


Hag =a(m"), or pa = a(m*)/q. (58) 


This is just the first scaling relation in (20). For the cases with m? = 5/4,0 and —3/4, 
we have peg ~ 2.784531829, 2.265193164 and 1.737772548, respectively. It is easy to see 
that they are very close to corresponding ones obtained by using the shooting method 
in section 3. This agreement also provides a piece of evidence to support the results we 
obtained by shooting method to solve the equations (11). 

It seems impossible to get the second scaling relation in (20) by an analytical method 
since it happens in the condensed phase with large condensate. In this case, the back 
reaction of the vector field can not be neglected, and we have to solve the equations (11) 
without any approximation. By numerical analysis, we find that when q is very large, the 
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In(u ,)=-0.33 LIn(q)+0.927 


Pe C=0.4791m7+ 0.0492m*-5.0*10~+ 


Inu.) 


-0.5 l 2 3 -0.353 07 06 05 04 03 02 01 0 
m m 

Figure 25: The fitting curve for the second relation in (20). In both figures, the red circle 
and blue line are the numerical results and fitting curves, respectively. The left plot shows 
the scaling relation between In(ue2) and In(q) in the case of m? = —3/4. It shows that 
there is a well-behaved linear dependence when In(q) is larger than 2. In the right figure, 
we plot ¢ with respect to different values of mass square. The numerical results are fitted 
well by a second order polynomial with a small constant term. 


critical chemical potential for the phase transition from the condensed phase to the normal 


phase has a power law form as ; 
ben x qg. (59) 


A typical case with m? = —3/4 is shown in the left plot of figure 25. The fitting of ¢ for 
different m? is shown in the right plot of figure 25. One can see that the relation between Ç 
and m? is fitted well by a second order polynomial with a very small constant term. This 
constant term in fact is less than our numerical error. Thus if neglect this constant term, 
we can obtain the following fitting relation 


C(m?) = 0.4791m? + 0.0492m*. (60) 
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